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Abstract . A general Hamiltonian for a rotating moon in the field of the 

earth is expanded in terms of parameters orienting the spin angular momentum 

relative to the principal axes of the moon and relative to coordinate axes 

fixed in the orbit plane. The effects of elastic distortion are included as 

modifications of the moment of inertia tensor, where the magnitude of the 

distortion is parameterized by the Love number k^. The principal periodic 

terms in the longitude of a point on the moon due to variations of the tide 
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caused by the earth are shown to have amplitudes between 379x10 and 176x10 
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with a period of an anomalistic month, 370x10 and 172x10 with a period of 

one-half an anomalistic month and 274x10 ^ and 976x10 4 with a period of 

one-half of a nodical month. The extremes in the amplitudes correspond to 

rigidities of 8xl0 + ^ cgs and 2xl0 + ^ cgs respectively, the former rigidity 

being comparable to that of the earth. Only the largest amplitude given above 

is comparable to that detectable by the projected precision of the laser 

ranging to the lunar retroreflectors , and this amplitude corresponds to an 

improbably low rigidity for the moon. A detailed derivation of the free 

wobble of the lunar spin axis about the axis of maximum moment of inertia is 

given, where it is sh own that elasticity can alter the period of the free 
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wobble of 75.3 years by only 3><10 ' to 10 of this* period. Also, the effect 

of elasticity on the period of free libration is completely negligible by 
many orders of magnitude. If the moon's rigidity is close to that of the 
earth there is no effect of elasticity on the rotation which can be measured 
with the laser ranging, and therefore no elastic properties of the moon can 
be determined from variations in the rotation. 
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ii. . Tn-t ro'duc tio n 

Order of magnitude estimates of the periodic variations in the rota- 
tion rate of the moon due to the changing magnitude of the tide induced by 
the earth indicate marginal detectability of these variations by future 
laser ranging to the retroref lectors on the lunar surface. It is thus 
appropriate to determine more accurately the magnitude of the effects of 
elasticity on lunar rotation to see whether any must be eventually included 
in the reduction of the laser ranging data. If such effects are in fact 
measurable, perhaps a more important result will be the determination of the 
effective elastic properties of the entire moon with possible implications 
about the, nature of the interior. 

Three possible effects of elasticity are considered here. The first 
is the above mentioned variation in the rotation rate due to tidal changes in 
the inertia tensor. In Section 2 the largest term of this variation is 
calculated directly from the tidal variation in the moment of inertia about 
the spin axis. This serves as a check on the general theory described in 
Section 3 from which all the perturbations of the lunar rotation, including 
those due to elasticity, can be found to arbitrary order. The amplitudes 
of the three largest terms of the angular displacement of a point on the 
lunar surface from its mean sidereal position (corresponding to uniform rota- 
tion) are determined to Section 4 for effective rigidities corresponding to 
those of steel and aluminum. 

For the second effect of elasticity to be considered, the theory devel- 
oped in Section 3 is used in Section 5 to determine the period of the free 
wobble as modified by the presence of the earth and to demonstrate that the 
effect of elasticity on the free wobble is negligible. This could be antici- 
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pated beforehand from the functional, dependence of the modification of the 
wobble period of the Love number k 2 (see Munk and MacDonald, 1960, for a 
discussion of Love numbers) and the spin rate ip and a comparison with the 
modification of the earth's Chandler wobble period. 

The period of the free libration in longitude is considered in Section 6, 
where any alteration of the period by elasticity is quickly dismissed. Even 
a small change in the free libration period is important because of the exis- 
tence of forced librations whose periods are very close to the free libra- 
tion period (Jeffreys, 1957; Eckhardt, 1970; Williams et al . 1973). A few 
percent change in the free period can lead to a change in the amplitude of 
the near resonant forced libration of a factor 2 or 3. But no measurable 
alteration, even with the amplification, is evident. 

Section 7 is a summary of the results where it is pointed out that 
only the largest term in the variation of the rotation rate is possibly 
of measurable amplitude and that only if the moon has what is perhaps an 
unusually low rigidity. 

2. Tid ally Induced Variations in . t he S p in __Rate 

If ip is the angle between the axis of minimum moment of inertia and 
the ascending node of the lunar equator on the orbit plane then with suf- 
ficient accuracy for the short term effects considered here, 

d 2 ip '^dijidC m 

, 2 ~ C dt dt U ' 

at 

for a tidal variation in the moment of inertia C about the spin axis (spin 
angular momentum is conserved). The time variation in C is determined from 
the variation in the tidal mass distribution. 

The magnitude of the tide at a point R on the lunar surface is given by 
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AR = h 0 —S' a- P 9 (cos0")' (2) 

z 3 e z 
r g 

where R is measured from the lunar center of mass, P 2 s t ^ ie Legendre poly- 
nomial, cos0" = R*r/(Rr), r being the position of the earth relative to the 
moon, G is the gravitational constant, g is the surface gravity on the moon, 
a is the lunar equatorial radius. M is the earth mass and h„ is the Love 

6 w Z 

number defined by (Munk and MacDonald, 1960) 

h 2 = — MV ‘ (3) 

2 p § a e 

for a homogeneous sphere, where p is the coefficient of rigidity and p is 
the lunar mean density. The tidal mass per unit area of the lunar surface 
Am/AA = pAR and 


C - C n + 77 a 2 sin 2 0dA 

0 AA e 


where 0 is the polar angle measured from the spin axis. 

Part of the mass in Am is determined by the fluid Love number (p = 0) 


and corresponds to the tide at the mean earth-moon separation. However, 
this constant tide vanishes in dC/dt, and h 2 defined by Equation (3) is 
appropriate for the varying tide. Then 


dC d Am 2 . 2 OJA 

77 = 7 “ 77 a sm 0dA 

dt . dt AA e 
A 


17 f 2 71 M a 6 

- 3h 2 P 2 (cos0") sin^0d0d<J> 

0^0 r 


where is the lunar mass, <f> is the azimuthal spherical polar coordinate 
of position on the moon measured from the axis of minimum moment of inertia 
on the lunar equator. Substitution of ma^p, dr/dt = nae sin f/ /l-e 2 

with n, a, e, f being respectively the lunar orbital mean motion, semimajor 
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axis , eccentricity and true anomally, cos0" = sin0cos<fi and performing the 


integration in Equation (5) yields 

5 


= -2 h M -SL. 
dt 5 2 « 3 

a 


e n sin nt 


( 6 ) 


4, 4 X 


where only the first term in the expansion of (a /r ) sin f in the mean 
anomaly (M = nt) has been kept, higher frequency terms being higher order 
in e. 

Substitution of Equation (6) into Equation (1) , with dtp/ dt and C replaced 

2 

by their mean values n and 0.4 M a g respectively, gives 


d 2 if> 3 2 a e M © . 

— 2 = ~~2 n — 3 ~ m - h 2 e sin nt 
dt a 


(7) 


which upon integration yields an amplitude of variation in ip from the 


mean (uniform rotation) value of 

3 


. a 

1 

2 3 

a 


M § |3V9xlO 

M 2 | 1V6x10“ 


«-3 
2 
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with a period of the anomalistic month. The two numerical values correspond 

2 

respectively to a lunar rigidity like that of steel (y = 8 x 10" dynes /cm ) 

2 

and aluminum (y = 2 x 10" dynes /cm ). The rigidity of the earth is slightly 
below that of steel (Munk and MacDonald, 1960), so one might expect that 
of the moon to fall within the above extremes. 

The single term in evaluated above is expected to be the largest 
effect of elasticity on the lunar rotation. Information about other terms 
in the rotational variations requires a more general theory, which is 
described below. This theory is sufficiently general to include rotational 
distortions and deviations from principal axis rotations and is used to 
investigate the free wobble in Section 5. 
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- 3. . General Theor y 

Parts of the following development were used in earlier investigations 
of the rotations of the moon. Mercury (Peale, 1969) and Venus (Goldreich and 
Peale, 1970). In those applications the Hamiltonian did not include the 
effects of non principal axis rotation which are necessary for an investiga- 
tion of the free wobble. These effects are added here, and an error in the 
definition of variables is also corrected here. 

A complete development of the variational equations and the form of the 
Hamiltonian is beyond the scope of this paper. Let it suffice then to out- 
line the procedure and then write the final form of the Hamiltonian and varia- 
tional equations. 

The origin of coordinates is at the center of mass of the rotating body, 
and the Hamiltonian is’ the sum of the rotational kinetic energy and the poten- 
tial energy due to external gravitational fields. The translational kinetic 
energy and central terms in the potential energy do not contain the coordi- 
nates and momenta associated with the spin and orientation and are therefore 
suppressed. Body fixed axes designated by the unit vectors i, j, k are the 
principal inertial axes. A second set of axes designated by (I, J, K) , are 
fixed in the orbit plane of the disturbing body with K being normal to the 
orbit plane. The earth-moon orbit precesses on the ecliptic plane so I is 
chosen to be along the mean ascending node on the ecliptic of the earth's or- 
bit relative to the moon. 

The (i> j> k) system can be oriented with respect to the (I, J, K) by the 

ordinary Euler angles (<J>,9,i|0. The angular velocities S , S , S can be 

x y z 

expressed in terms of the Euler angles (Goldstein, 1950) and generalized momen- 
ta p^, pg, p^ conjugate to <j>, 8, ip can be used to express the kinetic energy 
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part of H with the resulting variational equations being canonical. 

The set of variables (<f> , 0, P^> Pq> P^) is not convenient for expres- 

sing the Hamiltonian or for interpreting the variations. The following series 
of transformations is therefore effected: 

(<f>» 8, V, Pq* P Q » P^) -*■ (a, n, K*a, J*a, k-a, j *a) 

(a, n, K*a, ft, k*a, to) , 


where a = aa is the total spin angular momentum, a being a unit vector, n is 
the angle between Kxa and kxa and is thus a measure of the rotation about the 
spin axis. The variables ft and to are defined by 


cosft = 


-J*a 

/1-F>a 2 


I* a 

sinft = -rrrzzz: 
/l-K-a 2 


-j*a j *a 

cos = ■■ 1 ~ ■ sin = - 

/l-k* a 2 /l-k-a 2 

ft is thus the angle between Kxa and I and to is the angle between kxa and i. 
The above choice of variables conveniently orients the angular momentum 
relative to space and body axes . and describes the rotation. The variational 
equations in terms of the final set of variables are 


da = __3H 
dt 3n 


dn _ 

3H 


K-a 

3H 

k-a 

dt 

3a 


a 

3K* a 

a 

dK-a 


1 

3H 

K- a 

3H. 

dt 


a 

3ft 

a 

3n 


dft = 1 3H 
dt a 3K-a 

<%-a i 3 h 4 fe-a aH 

dt a 3w a 3q 


( 10 ) 
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d(j _ 1 3H 
dt a 3k- a 


( 10 ) 


The development of the potential part of the Hamiltonian follows that of 
Kaula (1961). The terms in the expansion used by Kaula are each expanded 
to second order in /l-k*a^ such that the variational equations are correct 
at least to first order in this quantity. This is sufficient accuracy, since 
all sizable solar system bodies are expected to be rotating nearly about a 
principal axis (k*a ~ 1) , driven there by energy dissipation. The develop- 
ment of the kinetic energy part of the Hamiltonian follows that of Peale 
(1969) except now the allowed distortion of the elastic body introduces incre- 
ments in the components of the inertia tensor which are determined by rota- 
tional and tidal distortions. As such, these distortions will depend on the 
variables used in Equations (10) , but this functional dependence of the distor 
tions is considered only after the differentiations have been performed on the 
right hand sides of Equations (10). The Hamiltonian is 


H = | a ' [I] - ya + V (11) 

where [I] ^ is the inverse of the inertia tensor, jj is the precessional angu- 
lar velocity of the orbit plane and V is the potential. 


a = ot[(i*a)i + (j*a)j + (k*a)k] (12) 

V. 

ja = -y[siniJ + cos iK] (13) 

where i is the inclination of the lunar orbit to the ecliptic. 

Substitution of Equation (12) and (13) into Equation (11) , use of Equa- 
tions (9) and expansion of the terms in V gives the final general form of H 
expressed in the variables of Equation (10). 
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a 

( (l-k-a 2 ) 

r , „ 2 • 2 q , 2 

2t 

~Wl 

1 2 

1 22 1 3 3 1 2 3 + I 11 I 33 ~ I 13 J + 

I ll I 22 " 1 12 j 


U-fc’a > r 2 2- 

2 C ° s2a) I H I 33 ~ I 13 ~ I 22 I 33 + I 23 


+ (1-k-a 2 ) sin2co[l 23 I 13 - I 12 I 33 


a 


- 2/l-k-a 2 sin™ [l 12 I 23 - I 22 I 13 ] 


-.‘2/l-k-a 2 cos™ [l 13 I 12 - I U I 23 ] 

+ ya^-sini /l-K*a 2 cos 0 + cos i K>aJ 

00 /a £ £ °? 

1 (-“) I l l G (e) . 

£=2 ' ' m=0 p=0 q=-=» £pq 


-®vs 


m=0 p=0 q=-“ 



COS V. 

£mpq 

£-m even / . \ 

( sxn v 
1 £mpql 

c 


+ S ) l 

£m 

. 

£mj 


sin v. 

£rapq 

/—cos V. 

t-n odd 1 ‘“H 


£.— m even ■ 


£-m odd - 


+ /l-k*a 2 A. (K*a) 
“ ~ £mp “ “ 


£m 


+ /l-k* a 2 B (K*a) 

~ ~ xmp ~ 


Zm 


COS<V impq +, ’ > 


Sin<V <l I n P q +n> 


COs(v impq- ,l) 


sin(v -n) 
£mpq 


£-m even 


sin(v impq +r,) 


£-m even-. 


+ S 


£m 


£-m odd 
\ £-m even 


l-cos(v„ +n) 

£mp q ‘ £-m odd J 


+ S 


£m 


£-m odd 


1 sin(v 0 -n) 
£mpq 


-cos(v. -n) 
£mpq 


£-m evenn 


£-m odd 


+ (l-k-a 2 )C. (K-a) 

~ £mp ~ 



C0S<V impq +2,l) 

£-m even 

Sln<V impq +2 " ) 

£-m even - 

c 


+ S 



£m 


£m 




sin(v tmp, +2n) 

^£-m odd 

- ct>s(v lmpq +2r,) 

£-m odd J 
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+ U-k-a 2 ) %,<K-a> 

l COS(v £mpq- 2ri) 
C lm ) 

£-m even 
; + S £m 

1 

sin(v. -2n) 

£mpq 

£-m even- 

. 

■ 


_ ( Sln(v £»pq- 2n) 

L-m odd 

cos(v. -2n) 

£mpq 

£-m odd - 


(14) con't. 

where I., are the components of the inertia tensor, |l| is the determinant 


of the inertia tensor, G, M^, M^, a, a g e have been defined previously, C^, 
are the ordinary harmonic coefficients in the expansion of the lunar 
gravitational field defined by 


(£-m):(2-6£ m ) 

’ 

r 

M a £ ^ ( £+m) I 

< 



(cosing 

P (r' ,6' ,<J>')r ,Jt+ t’^Ccose')*/ [ 

(sinm<f>) 


sin<j>' d9' d4>'dr' (15) 


the integration being over the entire volume of the moon. G. (e) are expan- 

£pq 

sions in the eccentricity and F^ (K»a) are functions which are both defined 

and tabulated by Kaula (1966) with (K*a) replacing cos i . A ^ m p» B £ m p» C £ m p ’ 

D. , E. are defined in terms of the F„ , and are given in Table I for 
Amp £mp Jimp 

£ - 2. Finally 


v. = ( £— 2p ) (r-fi+y) + (£-2p+q)M 4- m(iT-Ti-u) 
£mpq 


is equivalent to Kaula’s definition (1964) of this angle, but the variables 
except for the mean anomaly M, are not those used in that work. We have 
used the orbit plane as a reference whereas Kaula has used the instantaneous 
equator plane. In addition, non-rincipal axis rotation forces us to use 
the vector a parallel to the spin angular momentum for body orientation 
rather than the principal axis coincident with the spin vector. This intro- 
duces n and to into v„ , where in the limit of principal axis rotation 

£mpq 

(k*a = 1) n + to = ip locates the axis of minimum moment of inertia from the 
node of the equator on the orbit plane. The angle y locates the pericenter of 
the earth's apparent orbit about the moon relative to the X axis, which is 


- 12 - 


along the orbit node on the' ecliptic. 

In the moment of inertia tensor, the products of inertia come from 
tidal and rotational distortion of an elastic moon and are therefore small 
compared to the diagonal terms. The diagonal terms have small increments 
as well due to elastic distortions. This suggests that we write the diagonal 
terms 

I-. = A + AA 

I 22 = B + AB (16) 

I 33" C+AC 

where A, B, C are the principal moments of inertia in the order of increas- 
ing magnitude and AA, AB, AC are the contributions from elastic distortion 
of the same order as the products of inertia. The determinant |l| can then 
be expanded in the equations (after partial derivatives of H are taken) and 
only first order terms in AI../I.. and I.. /I., kept. This simplifies the 
equations of motion considerably and sufficient accuracy is ensured by the 
smallness of the increments. It will also be noticed that elastic distortion 
introduces a time varying increment into the harmonic coefficients and 

S^. However, since these coefficients are themselves very small and are 
preceded by another small coefficient, the effect of the small elastic varia- 
tions in and are second order or higher in small quantities and can 
be ignored. Elastic distortion has its major influence in the kinetic energy 
terms except perhaps when the secular effects of dissipation are considered. 

The increments in the moment of inertia tensor are evaluated by compar- 
ing the potential from the elastic redistribution of mass with the terms in 
the expansion of the moon's gravitational field depending on the second moments 
of the mass distribution. Both the rotational and tidal distortions are 
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caused by harmonic disturbing potentials of the second' d'egree (e.g., see 
Munk and MacDonald, 1960). For rotation, the disturbing potential per unit 
mass is 


v: - RS 


2,2 2 2 

— P 2 (cos6’) - 


(17) 


where S is the spin angular velocity, R is the position of the field point 
relative to the lunar center of mass and cos 6 ’ = R*S/(RS). The central 
term in Equation (17) can be absorbed by the central term of the general 
lunar field and will not be considered further. The tide raising potential 
is given by 


- GM„T 


V’ = 
T 


3 - P 2 (cos0") 


(18) 


where cos 6 " = R*r/(Rr). Love (1944) has shown that the increment in the poten- 
tial at the surface of a spherical body distorted by a spherical harmonic 
potential is proportional to that distorting potential and falls off exterior 
to the body as r where n is the degree of the disturbing potential. 

The external potentials due to the lunar mass redistributed by rotation and 
tides are thus 


k 2 a 

V R = T 


a 5 / „ 

t(i» 


s ) 2 - ± r 2 s 2 


) 


GM a o 

V = -k — — — 
V T *2 3 5 

r R 


3 (R-r)' 

2 272" 

r R 


(19) 


( 20 ) 


where k 2 is the Love number defined by 


k 2 


3/2 


1 + 


19 _u_ 
2 Pga 


( 21 ) 


for a homogeneous sphere. If R, S, r are written in terms of their components 
in Equations (19) and (20) and compared with a similar development of the 
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second order term in the' lunar gravit'afeionai- field, the two representations 
of the field are equivalent and the trace of the inertia tensor invariant 
before and after distortion only if the increments in the tensor components 
are given by 


AT 


k 2 ae 5 S 2 


ii 


3G 

5 


- i 7) 


k a k a 

f / e „ „ / e 

T ij 3G S i S j " 5 M ® r i r j 


i 

i * J 


( 22 ) 


We can separate out that part of the rotational increment which appropriate 
to the fluid Love number by assuming principal axis rotation (S^ = = 0, 


“ S) . These we shall include in the permanent moments A, B, C and write 

i - 1,2 


ii 


33 


,5 ,5 

k a _ k a 

2 e _ s 2 _ 2 e M r 

3G 1 5 M ^ r -- 


~ k 2 a e 5 .2 2. ~ k 2 a e M 2 

3G (S 1 + S 2 J 5 M « r 3 


(-if?) 


(23) 


where AI„ (i j) remains unchanged. The tidal increments will also have a 
component appropriate to the fluid Love number. However, since these involve 
the coordinates of the earth, the tidal terms can be expanded in terms of the 
variables and functions used in the Hamiltonian and it is a simple matter to 
remove the constant terms from these expansions such that the appropriate k£ 
for the remaining terms is that defined by Equation (21). The above develop- 
ment can be applied in a straightforward manner to the effects of elasticity 
on the lunar rotation. 


4^ AEEii£§tion_of _General_Theory_to_Fluctuations_in_the_Rotation_Rate 

For our purposes here, it is sufficient to assume principal axis rotation 
and to ignore the effects of the precessing coordinate system, since these can 
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be added at any time. Principal axis rotation eliminates the rotational 
increments in the inertia tensor, and since we have neglected second order 
effects, we need consider only the kinetic energy term of the Hamiltonian. 
Both n and m are undefined when k = a, but their sum locates the axis of 
minimum moment of inertia from the node. Hence, we wish to determine 


4r = iim — (n+w) = lim + “ (Hs*a) 

dt dt (k-a*l)L 3a a 


m 1 

3(k‘§) J 


(24) 


k* §->1 

where 8H/3(K*a) = 0 since only the kinetic energy terms are involved. The 
limit is taken after the differentiations of H. From Equations (14 and (24) 


ii - “ Tt T t 2 i 

dt JT[ [I 11 I 22 ~ X 12 1 


(25) 


Using Equation (16) and keeping only the first order terms in the expansion 
of both the numerator and denominator, we find 


djp _ a L Ac\ 

dt c\ c / 


(26) 


where AC is here the tidal increment. The angular part of AC(tide) is the 
same as that in the coefficient of C£q in the general potential of the moon. 
Hence, we can lift the coefficient of C£q from Equation (14) and write 
. . / o k„M a 5 2 

i? •! 1 + f-r^ V lG W e V a) • 

\ Ca p=0 q 


cos 




(27) 


2p) (H-fi+y) + (2-2p+q)M - e 


20pq 


where e 20 pq * s a Ph ase lag due to dissipation (Kaula, 1964) and is included 
here only for completeness. If we combine terms of the same frequency and 
use Cassini's laws to set ft = r, we have 
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dj» _ _a 2 _a_ k 2. M ©' 5 v 

dt C 3 „2 3 a e 

L> 3 


| 2F 200 G 20-l cos [ 2y+M] 
+ 2F 200 G 200 cos [ 2y+2M] 
+ 2F 200 G 201 COS ^ 2y+3M 1 
+ 2F 200 G 202 cos t 2 Y+^] 

+ 2F 201 G 212 COs[2M] 

+ 2F 20l G 211 COS ^ I 


(28) 


The last two terms depend on orbital position relative to the perigee and 

are hence related to the changing magnitude of the tide. The first four 

terms depend on position of the moon relative to the node and are related 

? 

to the latitude of the tidal maximum. With I^qq = - 3(1-K*a )/8, F 201 = 


5 2 


9 2 


4 4 K ’ a * G 200 1 ~ 2 6 ’ °212 4 6 ’ G 211 = 2 6 ’ = cos ( 6 °41')» 

e = 0.0549, integration of Equation (28) yields the following amplitudes 

in A for the three largest terms: 


Amplitude 
379x10 
176x10' 

370x10' 
172x10' 

2 74x10" 4 
976xl0 -4 




Period 


Anomalistic month 


half of anomalistic month 


nodical month 


The values in each pair correspond respectively to lunar rigidities like that 
of steel and aluminum. The term with the period of the anomalistic month is 
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just that determined in Section 2, and* its amplitude agrees with the one 
determined there. 


5. Free Wobble 


The wobble of the lunar spin axis about the axis of maximum moment is 
determined by the last two of Equations (10), In addition to the contribu- 
tions to d(k*a)/dt and dco/dt by the kinetic energy terms, those terms from 
the potential part of H which are constant or depend only on to must also be 
retained. All other terms in the variation of to and k*a will have phases 
and amplitudes determined by the forcing term and therefore do not contribute 
to the "free" wobble. Let us first select the terms which must be retained 
in the potential part of H. 

From Equation (14) a general term has arguments 


V 

impq 

v 0 + n 

impq - 

V. + 2n 

S-mpq - 


(29) 


But from Cassini’s laws fi = rr and n + w»^ = M+ ^ where ip is the 

o o 

value of ip at perigee. For completely damped librations, \p° = y completes 
the law of stable synchronous rotation. These conditions can be applied only 
after the partial differentiation of H in the equations of motion, but can 
be used beforehand to select those terms which will depend only on to. Elimi- 
nating n, and from the arguments by Cassini's laws leaves 


(£-2p-m)y + (£-2p+q-m)M + mr 
(£-2p-ntfl)y + (£-2p+q-m+l)M + mu + m 
(£-2p-m+2)y + (£-2p+q-nH-2)M + mir + 2to 
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corresponding to the five arguments in (29). The only permitted arguments 
for 1 = 2 which are constant or depend only on <o are 


V 2010 


constant 

v 2200 


constant 

V 2110 

+ n — * 

-0) 

v 2100 

- n — *■ 

+10 

V 2210 

+ 2n -»■ 

-2io 

v 2020 

+ 2n -*■ 

-2to 

V 2000 

- 2n 

+2(o 


Selecting only these arguments from the Hamiltonian, the variations of 
(o and k*a assume the forms 


d\p 

a) 4 a\ *** 

(i-^) 

2k* a / 
_L Z. f 1 

— M 

dt 

2 | A y 1 ~ A J B 

l 1 ‘) 

C 

C /\ 

+ 

2GM M 9 

3 \ *-2 (E 201 C 210 C 

aa 

4- F 

20 220 

G 200 C 22 ) 



a(l-2k*a 2 ) ''"23 , a 

- “ — — COSO) + 

(l-2k* a 2 ) 

I 13 . 

• r sinco 


+ 

/ o BC 

v ±-k’ a 2 

/l-k* a 2 

AC 



1 


2GM„M, 


P 2. 

— a k* a 


aa 


3 0 

e “ 


(30) 


C 221 G 210 C 22 + 


+ C 202 G 220 C 20 + D 200 G 200 C 2ol ^ cos2a) + ^-a 4r sin2co 


h 

~ AB 


d -’ 3 . /— 2 P 13 123 • 1 

, . = -ak*a vl-k’a^ I— — cosco — smu 

dt - - (_ AC BC J 


f i , 2. 12 . 

+ a (l-k* a ) -TT- cos2u 

** ** a K 


) a , 2. 

if 

ab\ 

if 

aa\ 

j j (l-k* a ) 

K 1 - 

V 

■ *( ' 

-)J 


( 31 ) 



-19- 


2GM M 

2 ,, , 2 n r 

3 a e ) < 

a a L 


C 221 G 210"22 


C "" + ( C 202 G 220 


+ D 200 G 200 


)C 2Q k*gj|sin2w 
(31) cont'd. 


In Equations (30) and (31) the terms involving C^, and S 22 have been 
omitted since they involve products of inertia which are only induced by 
elastic deformation and are preceded by a small factor. Also the AI„ 
will involve only the rotational deformation since the tidal increments 
involve terms which contain to only in the sum q + to which is transformed 
to M + i|) by Cassini's laws. Hence, there is no way to isolate a term depend- 
ing only on to in the AI^ due to tides and the tidal distortion terms in the 
wobble classify as forced terms. For the rotational increments it is 
sufficient for first order accuracy to write 


S i ~ I 


(32) 


and to set I=A=B=C=0.4Ma e in the denominators of Equations (32). 
Substitution of the above form for into Equations (22) (minus the tidal 
contributions) and use of Equations (9) gives 

AA \ 2 , 2, . 2 

— = C(i*a) = C (1-k* a )sin oj 

AB r /4 \ 2 . 2. 2 

— - £ (j • a) - £(l-k*a )cos u> 

AC 2 . . 2. . 2. 

— = ? (i* a +j • a ) = -;(l-k*a ) 


— ■ = ?(i*a)(j*a) = C (1-k* a^) sin2m/2 


( 33 ) 


13 


= £(i*a)(k*a) = -?/l- k*a2 ,k»a si 


m« 


^23 ' 

— = £ ( j • a) (k* a) = -?/l-k*a 2 k* a costo 


where 
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5 *2e 


6 GM. 


(34) 


with a = 1^ being used in Equation (34). 

Using the expressions in Equations (33) in Equations (30) and (31) we 
arrive at the expressions 


dcj 

dt 


= J^K 1 + K 2 cos2tijJk* a 


d(k*a) 0 

dt = -K_ ' (l-k*§ )sin2to 


(35) 


where 


h--k[ 


B-A _ 2 (C-A) 
B C 


+ 2Q 


] 


2G V 


aa 


3 ^ a e 2 [" E 201 G 210 C 20 + E 220 G 200 C 22 


K a < B ~ A ) , 2GM A , 2f c f c , fc 

K 2 ~ 2 AB + 3 e [ C 221 G 210 C 22 (C 

aa *- 


(36) 

202 G 220 + D 200 G 200 ) C 2o] 


2GM M 

K 2 = 2 ( AB ) + 3^ a e [ C 221 G 210 C 22 + (C 202 G 220 + D 200 G 200 ) C 20 (k ‘ a) ] 

aa L J 

In equations (36) , c/A = ?/B = £/C was assumed for first order accuracy. 

The expressions for K 2 and K 2 differ only by the factor k*a multiplying some 
of the terms. However, d(k*a)/dt contains l-k*a as a factor so it is consis- 
tent with the first order accuracy to set k*a = 1 in in which case K 2 = 

K r 

Equations (35) are more easily solved in terms of the variables (i*a), 
(j*a) which are related to k*a and u) by Equations (9). This transformation 
(with K 2 = Kp yields 

d( ~‘ §) = (K +K 9 )(j*a)(k-a) 
dt 1 1 ~ 
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d ( j * a) 

— Jl (K 1 -K 2 )(i-a)(k*a) 

= 2 K 2 (i*§)(j-a) 


(37) 
cont' d. 


which are just a form of Euler’s equations for force free rigid body motion. 
For small deviations from principal axis rotation (i*a, j«a << 1, k*a « 1), 
the solution is immediate with 


i*a = ^cosCAk^) (K 1 -K 2 )t + cf^) 

/K -K 

i-s * Jvqzq + ♦]> 


(38) 


giving the expected elliptical motion of a about k. The shape of the ellipse 
and the wobble frequency follow from the values of and k 2 which are obtained 
with 


'202 

ii 
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co|t— • 
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220 

= -^(l+K- S ) 2 

221 

= - 3 
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212 

3 2 

= -|(l-3K-g ) 

'200 
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220 

= -^(1+K-a) 2 

201 
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N3|bJ 

F 

201 

= -|(l-3K*a 2 ) 

'220 

1 

"2 

F 

220 

= -|(l+K-a) 2 


(39) 


where K*a = cos(6°41'). All of the G functions in and K 2 differ from 1 only 

2 

by terms of order e so will not be listed explicitly. Placing numerical 

values into Equations (39) and expressing moment of inertia differences in 

C_b 

terms of a', 3', y' (correct to first order) where a’ = -——,3' = (C-A)/B, 

A. 

y’ = (B-A)/C, we find 
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K l +K 2 ~ ^y'-g'+cj. - . O01 *|j ja'+S , '+Y '} 
K x -K 2 ~ ^|-g+£ | - 1. 48i)j|a ,: fg ’+y 


(AO) 


where the second terms on the right hand sides of Equations (40) explicitly 
demonstrate the influence of the earth (with n = ifi) . The effect of elasti- 
city is contained in the parameter £. That elasticity has essentially no 
effect on the wobble is seen by comparing the value of c, with a’, 6', y ’ . 
With 


a’ = 3.97x10 H 
6' = 6.27xl0~ 4 
y f = 2.30xl0 -4 
? = 1.5xl0 -7 to 6xio 


(41) 


we see that 4 influences the motion by only a few parts in 10 to a part in 
3 

10 and is likely to be comparable to or smaller than some of the neglected 
contributions to the free wobble from % = 3 terms in the potential. 

This lack of influence of elasticity could have been anticipated by the 
functional dependence of C and a knowledge of the change in the period of 
the Earth's Chandler wobble by elasticity. The ratio of £/g' is the important 
parameter as seen from Equations (40) . 


_L ~ 
6 ' ~ 


k 

B 



' 2 


n 


(42) 


o 

where n is the orbital angular velocity of a satellite near the surface of 
o 

the moon (or earth). The ratio in Equation (42) is about 0.35 for the earth 
with 4>In o = 1.5/24 and k 2 = 0.3, and one observes about this fractional in- 
crease in the period of the Chandler wobble over that for a rigid earth. For 

-4 -3 

the moon this ratio is seen from Equations (41) to be 2.5x10 to l x 10 
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The parameter C depends on the Love' number arid' on the square of the rotation 
period. The former quantity drops for the moon by an order of magnitude from 
that of the earth and the latter drops by three orders of magnitude. These 
two effects combine to reduce the influence of elasticity on the lunar wobble 
by a comparable factor from the influence on the Chandler wobble. 

Finally with K^+K^ = -3.98x10 i and K^-K^ = "2.48x10 %, the ratio of 
the major and minor axes of the elliptical path of a relative to k is 2.57 
with the long axis perpendicular to the axis of minimum moment of inertia. 

The wobble period is 75.3 years, which had been obtained earlier (Sekeguchi, 
1970). This period would be increased by a factor of 2 if the earth were 
removed. 

6. Free Lib rat ion 


The existence of a forced libration whose frequency is very near the 
3 year period of the free libration motivates a check on possible alternations 
of the free libration period. However, we can quickly dismiss the effects 
of elasticity on this period. The major effect will be the tidal distortion, 
but the variations in the rotation rate discussed in Sections 2 and 4 are 
high frequency and will not disturb the libration. We are thus left only 
with the tidal torque arising from a dissipation caused phase lag. In the 
limit of small tides the net torque on a librating moon is just the sum of 
that on the permanent lunar bulge and the tidal torque. The latter is given 
approximately by 


T = -3k 


GM ® 2 a c 5 
2 2a 6 « 


(43) 


where Q is the specific dissipation function (MacDonald, 1964), The torque 
on the permanent deformation is determined from the £mpq = 2200 term in the 
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Hamiltonian with 2(ip -y) = 26 being the lib ration angle. With \p = , 


we have Cip = C(^ o ) = c (^ 0 _ y) since y = 0, or 


C 6 = 


~ G Hi G F (B-A) fi _ 3 k 2% . 5 

3 zOO F 220 C 6 3 2a 6 Q e 


If we make Q inversely proportional to frequency, then 
3 a e GM # G 200 F 220 (b _ a) 

4 + 2 — ~rr + ~7 ~ 4 ■ 0 

Ca o o a 


(45) 


where Q q is the value of Q at a reference angular velocity <5 q . With the 

. 2 ' 
coefficients of 6 and 6 being oi and 2 E, respectively, the libration frequency 

o 


is just /to 2_^2 = u [l-^^/(2o) )]. For F 0 __ ~ 3 and G„_. ~ 1, we have that 
O O O ZZU Z UU 



01 


(46) 


o 

does not affect the period. A value of Q = 100 with 6 corresponding to a 

o o 

3 year period was used in evaluating £. That the only effect of tides is 
the ordinary torque due to a phase lag of the lunar response is verified by 
evaluating da/dt from Equations (10) and selecting those terms in AI_ with 
arguments which are integer multiples of (ip Q -y). All terms with these argu- 
ments cancel exactly to first order except for the effects of phase lags con- 
sidered above. 


7 . Discussion 

The results of the previous sections imply that the only possibly 
important perturbations of the lunar rotation due to elasticity are the perio- 
dic fluctuations in the spin angular velocity. Even these are so small that 
their measurement cannot be expected in the forseeable future. The ultimate 
range accuracy of the laser radar to the moon now anticipated -is on the order 
of a few centimeters. (P. Bender private communication, 1972) which is 
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comparable to the linear displacement on the lunar equator corresponding to 
0701 shift in longitude. The results of Section 4- show that the amplitude 
of the longitude variation from the mean is comparable to this value only 
for the low rigidity limit. Since siesmic wave velocities on the moon are 
comparable to those in the earth's upper mantle (Toksoz , . et al . , 1972) and 
the density of the mantle and moon are comparable, one infers comparable 
rigidities of the moon and earth. Since the smaller amplitude of Equation 
(8) is thus more likely appropriate to the moon, measurement of the effects 
of elasticity on rotation must await the development of the next generation 
of instrumentation perhaps requiring placement on the moon itself. 

On the brighter side these results indicate that elasticity can most 
probably be safely neglected in the reduction of the laser ranging data. 
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